Interactions between colloids induced by a soft cross-linked polymer substrate 
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Using video-microscopy imaging we demonstrate the existence of a short-ranged equilibrium at- 
traction between heavy silica colloids diffusing on soft surfaces of cross-linked polymer gels. The 
inter-colloid potential can be tuned by changing the gel stiffness or by coating the colloids with a 
polymer layer. On sufficiently soft substrates, the interaction induced by the polymer matrix leads 
to large-scale colloidal aggregation. We correlate the in-plane interaction with a colloid-surface 
attraction. 

PACS numbers: 82.70.Dd, 82.35.Lr, 68.08.-p 



Knowledge of the interactions between colloids con- 
fined to surfaces [TJ— [3[ or interfaces [4] is key to the 
fundamental understanding of many physical phe- 
nomena. For instance, interaction between colloids 
on surfaces can induce the self-assembly of ordered 
phases [5[ that may find applications in the engineering 
of photonic crystals [5-7]. Several aspects of cellular 
morphology |§[ Q , mechanical properti es p| [To| , mo- 
bility n, 111 and cell differentiation 9, 1211 have been 
found to be sensitive to the elastic properties of the 
environment. Experiments in this field often employ 
cross-linked polymer-gel surfaces [H, |To| 5 [ll[ as a model 
for biological tissues. We demonstrate that soft polymer 
substrates may actively induce non-negligible interac- 
tions between inanimate, micrometer-sized, objects. 
Analogous phenomena may affect as well the behavior of 
living systems of the same dimension and must be taken 
into account when the results of above cited experiments 
are analyzed. In this Letter we explore the nature of 
such substrate-induced inter-colloidal forces and show 
how they are modified by changing physical properties 
of the soft polymer gels. 

The experiments were performed using plain sil- 
ica colloids with diameter a = 1.16 ± 0.05 jam 
(Microparticles GmbH, Berlin) and a nominal 
density of 1.5-2 g cm -3 suspended in 100 mM 
tris(hydroxymethyl)aminomethane (TRIS) buffer at 
pH 8. Under these conditions the silica colloids have 
a net negative charge with a (^-potential of -42±1 mV. 
As soft surfaces we used 150 /am thick polyacrylamide 
(PAA) cross-linked gel films deposited on top of a 
microscope coverslip. The polymerization of the PAA 
was triggered by adding tetramethylethylenediamine 
and ammonium persulfate to a solution containing acry- 
lamide monomers and the cross-linker bis-acrylamide 
(bis-AA) in phosphate buffered saline buffer (PBS, 
pH 7.4). For control experiments we coated the same 
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silica colloids with positively charged poly-L-lisine- 
poly(ethylene glycol) (PLL-PEG, Surface Soultions, 
Switzerland) that is easily adsorbed onto silica and 
provides steric stabilization. PLL-PEG coated beads 
were almost neutral with a ("-potential of -2. 2 ±0.2 mV. 
Non-adhesive hard surfaces were obtained by coating 
the glass bottom of 8 mm diameter incubation wells 
(Sensoplate, Greiner bio-one) with PLL-PEG. 

Using a PAA concentration of 5 % we tuned the stiffness 
of the material by changing the bis-AA fraction. The 
observations of weak attractive interactions between 
colloids were qualitatively and quantitatively similar for 
gels with elastic shear modulus G' [l3[ between 55 and 
522 Pa, thus, we focus on a substrate prepared with 
4 x 10" 4 v/v bis-AA resulting in a G = 240 ± 5 Pa - 
we refer to this gel as the "soft" substrate. In contrast, 
colloids on "ultra-soft" gels obtained with bis-AA 
concentration of 1 x 10 -4 v/v and with G' = 12 ± 3 
Pa exhibited stronger interactions. On soft gels, sedi- 
mented colloids did not undergo large scale aggregation, 
however, we observed the formation of dimers, triplets, 
and small chains and clusters, as shown in Fig. QJa) 
and []Jb) f° r two different surface densities. On the 
ultra-soft gel, sedimented beads initially formed small 
aggregates that were still mobile and merged in larger 
ones. The aggregation stopped after about one hour 
as the reduced mobility of big aggregates kinetically 
hindered the process. A sequence of images of the 
aggregation is shown in Fig. DJc). The larger aggregates 
could be either amorphous or ordered in a hexagonal 
close-packed lattice [see Fig. DJd)]. The mobility of 
colloids within the aggregates was small and we only 
observed rearrangement of the beads located at the 
outer perimeter of the clusters. Using PLL-PEG coated 
colloids instead of plain ones, the aggregation still took 
place on the same timescale but the clusters showed 
liquid-like behavior with a continuous rearrangement of 
single colloids in the bulk of the aggregates. 

We measured the effective pair-potential of sedimented 
beads following the method proposed by Crocker and 




FIG. 1. Panels (a) and (b): microscope images of plain silica 
colloids sedimented on a PAA soft gel with different surface 
coverage. Panels (c)l-4: aggregation process of plain silica 
colloids sedimented on ultra-soft PAA gel surface at t = (1), 
t « 15 min (2), t « 30 min (3) and t w 60 min (4). Panel (d): 
highlight of an ordered aggregate on the ultra-soft gel. 



Grier [3] that is sometimes referred to as "blinking 
optical tweezers" (BT) [T^j. In the BT experiments, we 
positioned two isolated colloids at separations between 
1.5 and 10 /am using optical tweezers [16]. Then, just 
after releasing the colloids, we took movies of 15 s 
at 30 frames per second. The separation r between 
colloids was obtained by tracking their positions using 
conventional tracking algorithms [17j. After discretizing 
r with a mesh size of 25 nm we sampled the transition 
matrix elements that express the probability for the 
separation to evolve from the bin i to the bin j in the 
time between two consecutive frames. The stationary 
probability distribution of r is p\ = P%jP S j and the effec- 
tive pair-potential is given by V(r)/(k^T) = — log(p s ). 
The BT method has the virtue that it produces good 
statistics, but there are two known artifacts with this 
technique; optical interference 0, [H| and hydro dynamic 
interactions [3]. Superposition of colloid diffraction- 
images can result in errors in the determination of the 
inter-colloid distance r, and this can lead to systematic 
errors in the estimate of V(r) 0, [T^j. We corrected our 
data for these errors following the method proposed by 
Polin et al. jj. 
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FIG. 2. (color online) (a) Pair potentials of colloids sedi- 
mented on gels and glass surfaces measured with the blinking- 
tweezer (BT) or pair-distribution (PD) methods. Bare silica 
colloids sedimented on soft PAA gel surface measured with 
BT (A o) and with PD (Al solid line); PLL-PEG coated sil- 
ica colloids sedimented on a soft PAA gel surface measured 
with BT (B); PLL-PEG coated silica colloids sedimented on a 
ultra-soft PAA gel surface measured with BT (C); PLL-PEG 
coated silica colloids sedimented on a PLL-PEG coated rigid 
glass surface measured with BT (D, o) and with PD (Dl, 
dashed line) . The curves have been vertically shifted for clar- 
ity. Not-to-scale sketch of two beads partially included in the 
gel surface and interacting through excluded volume effects 
(b) and colloid-polymer adhesion (c). 



For the case of plain silica colloids supported by a soft 
PAA gel surface, the apparent V(r) shows two main 
features: a long-range attractive minimum with a depth 
of about 0.5 k^T and a sharp minimum for r w <r, 
as shown in Fig. Efa), curve A. Figure [2fa) curve B 
demonstrates that colloids covered with a PEG brush 
had a weakened short range attraction when compared 
to bare colloids on the same soft gel. In Fig. [2ja) 
curve C we show V(r) measured for PLL-PEG coated 
colloids on the ultra-soft gel [l9|. Here, the short-range 
attraction is stronger than for soft gels, while the 
long-range minimum is absent. In Fig. [2fa) curves A, B 
and C there is a maximum of the potential located at 
r « 1.2 a. This cannot be ascribed to screened Coulomb 
repulsion as the Debye screening length is estimated to 
be only a few nanometers [20]. As a control experiment 
we measured V(r) with BT for the case of PLL-PEG 
coated colloids sedimented on a PLL-PEG coated rigid 
glass surface. The potential, plotted in Fig. [2fa) curve 
D, does not show the short-range minimum observed 
on soft surfaces, although the long-range attraction is 
stronger. 

Because the long-range attraction in the effective pair- 
potentials for colloids was observed for both soft gels and 
rigid glass surfaces, it must have the same origin. The 
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form of the long-range attractive potential resembles 
the attraction observed between like-charged colloids 
in the presence of two confining walls [21]. As pointed 
out by Squires and Brenner the potential extracted 
using BT method arises from forces which may not have 
an underlying equilibrium potential. In particular, a 
hydrodynamic coupling between two moving colloids 
in the presence of a confining wall may produce an 
apparent effective attractive interaction similar to the 
long-range attraction that we measured. 
To separate this non-equilibrium hydrodynamic effect 
from the equilibrium potential between supported col- 
loids, we also evaluated V(r) from the radial distribution 
function g(r) obtained by measuring positions of colloids 
in equilibrium conditions using artifact corrected video 
microscopy [2j. For each sample we took a few hundred 
pictures at time intervals of 1 s of an area of about 
1.5xl0 -2 mm 2 with a surface coverage of 0.5-1 %. At 
these densities we can neglect many body effects and 
extract the pair-potential as V(r)/(k^T) « — log[g(r)]. 
We determined V(r) from measurements of g(r) for the 
cases of plain colloids settled on a soft gel surface and 
PLL-PEG coated colloids settled on a rigid glass surface. 
The potentials are shown in Fig. [2fa) curves Al and Dl 
respectively. In both cases, the long-range minimum is 
absent, consistent with the explanation that this attrac- 
tion was due to hydrodynamic coupling. The variations 
that we observe in measured potential energy between 
ultra-soft gels, soft gels and rigid surface correlates with 
differences in the colloid-surface interaction, as well as to 
a penetration of the hydrodynamic flow in the polymer 
network. 

The short-range attraction is present for the cases of the 
soft and ultra-soft gel surfaces, but is absent for the rigid 
glass surface as measured using both the non-equilibrium 
BT and equilibrium pair-distribution methods. These 
observations, taken together, indicate that the short- 
range colloid-colloid attraction is an equilibrium effect 
mediated by the polymer gel substrate. The same argu- 
ment applies to the repulsive maximum in the colloid 
pair-potential located at r « 1.2 a that is present when 
the colloids are on soft gels, but absent when colloids 
are on a hard glass interface, as shown in Fig. [2j We can 
deduce that the short-range attraction is responsible for 
the small aggregates on the soft gel substrates and for 
the large-scale clustering on the ultra-soft ones, where it 
became stronger. The weakening due to the PLL-PEG 
coating accounts for the liquid-like behavior of clusters 
on the ultra-soft gels. 

The range and shape of the short-range attractive 
well are compatible with a description in terms of 
depletion [22]. The PAA gel can in principle release 
free polymers in solution that may act as depletants. 
However, the absence of aggregation on bare glass in 
the same sample cell that contains gels with associating 
colloids excludes this possibility. Chen and Ma [23| 
predicted that a combination of excluded volume effects 
and polymer-colloids interfacial energy can give rise to 




-10 12 3 

z(jim) 



FIG. 3. (color online) (a) Experimental colloid-surface inter- 
action potentials for PLL-PEG coated silica colloids confined 
by PLL-PEG coated glass surface (A), plain (B, filled sym- 
bols) and PLL-PEG coated (Bl, empty symbols) silica col- 
loids confined by a soft PAA gel surface; plain (C filled sym- 
bols) and PLL-PEG coated ( CI empty symbols) silica colloids 
confined by an ultra- soft PAA gel surface. The solid line is 
a linear fit. Dashed and dot-dashed lines are quadratic fit. z 
is the direction normal to the surface. The curves are shifted 
vertically for clarity. In the insets: sketch of coupling between 
surface-colloid and colloid-colloid interactions for soft (b) and 
ultra-soft (c) gels. 



an interaction between colloids embedded in a polymer 
brush. A similar effect can occur in our system assuming 
a partial inclusion of the colloids in a weakly cross-linked 
polymer layer on top of the substrates. 
We explored this hypothesis by measuring the colloid- 
surface interaction energy. We tracked fluctuations of 
the distance z of sedimented beads from the surface by 
video-microscopy measurements [24[. The diffraction 
image of a colloid changes as a function of the distance 
between the bead and the focal plane of the objective 
lens. We took a set of a few hudred pictures of an immo- 
bile colloid by moving the focus with steps of 25 nm, then 
we associated an intensity I(z) to each reference image 
by integrating the pixel value within a circle concentric 
with the bead. Fitting the intensity as a function of z 
we obtained a calibration curve that is monotonic within 
an interval of a few microns. Working with the same 
illumination conditions we took movies of single freely 
diffusing colloids, then we measured I(t) for each frame 
and found the corresponding z(t) using the calibration 
curve. We sampled the probability distribution p(z) 
by taking a histogram of z(t). The colloid-surface 
potential is given by V(z)/(k^T) = —\og[p(z)]. To 
check the efficacy of the technique we measured V(z) 
for PLL-PEG coated colloids sedimented on PLL-PEG 
coated glass surface. As expected, the only potential 
attracting the bead towards the surface is the gravi- 
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tational one, as shown in Fig. [3fa) curve A. Fitting 
the linear region of V(z) we estimate the weight of 
the colloid mg ~ 2.6 fN, which is in good agreement 
with the nominal value > 3 fN. For plain silica colloids 
settled on the surface of a soft PA A gel, V(z) reveals an 
additional attraction, as shown in Fig. [3^a) curve B. We 
quantify the attraction with a parabolic fit that gives an 
effective spring constant of 2.1 x 10 -9 N m _1 . For the 
case of an ultra-soft gel, Fig. (3]^a) curve C shows that 
the colloid-surface attraction is even stronger with an 
effective spring constant of 3.6 x 10 -9 N m _1 . 
The attraction could lead to the partial inclusion of the 
beads into the gel surfac e 12 51 thereby accentuating a 
depletion-like interaction [23], as sketched in Fig. [2fb). 
Alternatively, a direct polymer bridging mediated by 
colloids-polymers surface adhesion can be responsible 
for the in-plane attraction, as shown in Fig. [2fc) (26| . 
Moreover, for the ultra-soft gels, the stronger colloid-gel 
attraction can result in a deeper penetration of the 
beads, which explains the stronger in-plane attraction 
either by increasing the excluded volume between two 
adjacent colloids or increasing the bead-gel contact 
surface. PLL-PEG coating on the colloids can reduce 
the polymer-colloid adhesion energy, thus the in-plane 
attraction, as found in V(r) measurements. The corre- 
lation between in-plane and colloid-gel forces is sketched 
in Fig. |3(b) and (c). The repulsion found for r « 1.2 a 
can be explained as the elastic response of polymer coils 



squeezed between two adjacent beads. We also measured 
V(z) on the soft and the ultra-soft gels using PLL-PEG 
coated beads, as shown in Fig. [3] curves Bl and CI. In 
the first case, steric stabilization does not affect V(z) 
while in the case of ultra-soft gels the potential is slightly 
less attractive. 

In summary, we demonstrated that silica colloids 
sedimented onto soft cross-linked polymer gel surfaces 
interact through a strong short-range non-hydrodynamic 
attractive potential that can be tuned either by changing 
the substrate stiffness or bead-surface properties, and, 
for the case of ultra-soft substrates, produces large 
scale aggregation, eventually leading to formation of 
ordered phases. The in-plane interaction correlates with 
a colloid-gel attraction that could result in a partial 
penetration of the colloids into the upper layers of 
the gels. All the experimental observations are consis- 
tent with the in-plane interactions arising from either 
depletion effects or direct polymer bridging. Further 
investigations are needed to determine which of these 
proposed mechanisms account for the aggregation of 
colloids on surfaces of soft polymer gels. 
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